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This is the final report on work performed under NASA Grant 
iio. NS6 1295. The grant is entitled "Dynamic Response and Input 
Identification". The work was performed at the University of 
Lowell, Lowell, Massachusetts. Research investigators were Dr. 
John A. McElman, Dr. John C. O' Callahan and Dr. G. Dudley Shepard. 
Graduate research assistants assigned to this grant were Chris- 
topher Gianaeopoulos, John R. Crischow and Chaur-Ming Chow. 

Dr. Robert W. Fralich was the grant monitor assigned by the NASA 
during the first phase of this study. This work is described 
under Part I of this report. Dr. Michael Card was the grant 
monitor assigned by the NASA during the second phase of this 
study. This work is described under Part II of this report. 
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1.0 Preface 


The research supported by NASA Grant NSG1295 had at the outset 
three major goals: 

1. To develop a general method for determining the response 
of a structure to combined base and acoustic random 
excitation. 

2. To develop parametric relationships to aid In the design 
of plates which are subjected to random force or random 
base excitation. 

3. To develop a method to Identify the Individual acoustic 
and base Inputs to a structure with only a limited number 
of measurement channels, when both types of excitation 
act simultaneously. 

A discrete method of analysis was chosen early In this study 
because it provided the flexibility necessary to analyze any type 
of structure that might be encountered. Study of both the response 
and the input identification problems led to the development of a 
computer program which requires as input only the modal parameters 
of the structure (natural frequencies, mode shapes, and damning 
ratios) and appropriate excitation or response spectra. This 
approach is ideally suited to be coupled with a finite element 
program. 

This study included both analytical and experimental work, 
and it is useful to emphasize the inherent contrast between these 
two approaches. In the theory-based analytical phase of this 
study, questions were first posed and then answered by the analyst. 
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whereas In the experimental phase, the questions were posed by the 
experimental data itself. The experimenter was then challenged to 
find answers to questions which were often unanticipated. The re- 
porting of this study has fallen naturally into two parts, one 
which describes the theory based approach and the other which des- 
cribes the experimental reality. 
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2.0 Symbol List 
English 


Greek 


a Fourier transformed input 

acceleration vector 

A cross section area of 

acoustic transmission line 

c speed of sound in iir 

C damping matrix 

Cx. coincident component of force 
and base acceleration cross 
power spectrum 

D spatial distribution matrix 

E elastic modulus 

f force vector 

g acceleration of gravity 

G power spectral density (PSD) 

H transfer function matrix 

I identity matrix 

L length of acoustic trans- 

mission line 

m mass; also modal degree of 

freedom index 

M mass matrix 

n physical degree of freedom 

index 

N(w) Fourier spectrum of noise 
source 

p pressure; also transfer 

function pole 

P(s) characteristic function for 
transfer function poles 

q modal displacement vector 

Q generalized force; also 

quality factor of resonance 

Q- quadrature component of force 
and base acceleration cross 
power spectrum 

S Laplace transform variable 

t time; also plate thickness 


8 frequency ratio ui/w i also 

plate aspect ratio ” 

2 

Y * coherence ratio 

r modal participation factor 

a(s) characteristic function 

C damping ratio 

n material damping loss factor 

V Poisson's ratio 
C(x,t}displacement function 

p mass density 

cross spectrum phase angle 
u frequency, radians/sec. 

(D^ natural frequency 

Subscripts 

a base acceleration input process 

f force input process 

M material 

p plate 

r response variable 

s specific; also shaker 

y plate deflection response process 

Superscripts 

* complex conjugate 

• time differentiation 

(~) (overbar) implies a modal vector 

or matrix 
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English (continued) 

U modal transformation matrix 

W(ui) Fourier spectrum of output 
signal 

X axial coordinate for acoustic 
transmission line 

X(w) Fourier spectrum of input 
signal 

y plate deflection response 
vector 

Y Fourier spectrum of response 

signal 

z transfer function zero 

Z(s) characteristic function for 
transfer function zeros 
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4.0 P;rt I: Theory-Based Models for Structural Response 

and Input Identification 

4.1 Introduction to Part I 

I 

The estimation of the vibration response of a continuous 
structure to a spatially distributed wide-band random excitation 
combines essentially two problems: the identification of the 
dynamic properties of the structure and the handling of the sta- 
tistics cf the multiple excitation and response variables. The 
inherent computational complexity of dealing with many variables 
which are both dynamically and statistically coupled often means 
that for real engineering structures the level of detail is 
limited by cost. 

At present solutions to the response problem use either the 

classic transfer function approach or the non-classical method 

of statistical energy analysis (SEA).^ Broadly speaking, classi- 
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cal methods which use continuous system models offer excellent 
spatial and frequency resolution but are limited to highly Idealized 
structures with simplified boundary conditions. Classical methods 
based on finite system models can handle more realistic engineering 
structures but the cost limit on the number of degrees of freedom 
in turn limits spatial resolution and frequency bandwidth. 

Finally, the SEA method offers wide-band frequency estimates at 
the expense of reduced spatial detail caused by the crude averaging 
of the spatial modes of response. 

This study uses the classical approach on a finite system 
model. It is assumed uiat the system is known in terms of its 
natural modes ond frequencies from either existing analytical solutions or 

b 


m 


from d finite element proyram. Dampiog ratio assumptions are 
made separately based on experimental data. An efficient compu- 
tational strateny called Pi'ogram flAND is then developed to 
perform both dynamic response anH input identification calcula- 
tions for a multi-degrec-of-free‘Jc(.( (MDOF) structure subjected 
to statistically coupled pt'essure and base excitation. 

4.2 Development of Program RAND 

4.2.1 Analysis 

The analysis proceeds in two steps. First, it is necessary 
to perform a signal analysis, i.e., to trace the ef-act of the 
infjut vector through the structure to the resulting output res- 
ponse vector. Next, a statistical analysis is made to determine 
the space and time correlations between the input and output 
vectors. 

Consider first the case of a continuous structure excited 
by a force process which is randomly distributed in space and 
time. The structure, is divided into n discrete elements. Ar 
n-dimensional force vector f(t) causes a displacement response 
described by an n-diniensional displacement vecior y(t). In the 
usual fashion tlic modal transformiation 

y = Uq (4-1) 

is used to trai'iSfont tiie orieinal oovertiino matrix difiereritial 
equation into tne uecoigjled modal equation 

I-lq + Cq -» Kq - <) (4-2) 

where U is the uiidauiped moOdl matrix, H, C. \ ar<* oirxjonal nwdal 
proper'ty metiroes, and itje t.odal rurcc veclot' f’ - U‘t. fuch a 
transfotiiiation existu wfienever C is proportional to or M; the 


transformation is approximately valid for non-proportional damping 
whenever the damping is light. 

For reasons of econoiny of computation it is desirable to deal 
with a reduced number of modal equations and modal coordinates 
equal to m, where m < n. Such a reduction is warranted, for ins- 
tance, whenever a mode does not significantly contribute energy 
to the ovt .1 response. So the matrices M, C, and K are of order 
(m X m), and U is of order (n x m). 

Since each mode is governed by a scalar equivalent to the 
matrix equation (4-2), a modal transfer function can be obtained by 
Fourier transformation for each modal coordinate. For Instance, 
for the p^” mode 


q = Q H (u)) 

^s,p 


(4-3) 


where p “ specific modal force. When all the mode 

response amplitudes are known the modal response vector q must be 
transformed back to the physical space vector y by the transforma- 
tion (4-1 ). 

All of the preceding signal operations that are required to 
relate the vectors f and y can be summarized in the signal vector 
block diagram shown in Figure 4-1 . 


1 n 

U H M U ' 




Fig. 4-1 Force Inoui; Transfer Ilatrix 
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This diagram uses the convention that the output vector is obtained 
by premultiplying the input vector by the transfer matrix. The 
matrix Fr(u) is a diagonal modal transfer function matrix consisting 
of the individual transfer functions Hp(io). The physical -space 
transfer matrix is seen to be related to the modal-space transfer 
matrix by the equation 

__-l T 

= U H M U' (4-4) 

Next the statistical relations which correspond to the signal 
vector block diagram can be determined by Fourier transforming 
related time domain correlation matrices. For instance, to the 
real coordinate transformation y = Uq there corresponds the fre- 
quency domain power spectrum relation 

Gy(a,) = UGp(u>) (4-5) 

and where the elements of G refer to the one-sided spectra cus- 
tomarily used in engineering calculations, and the single subscripts 
y and q refer to the y and q vector processes, respectively. This 
simple notation has been used rather than the more customary 
double subscript notation in order to avoid confusion later when 
more than one input vector process contribute to the overall res- 
ponse. For the complex transformation q = H(u)) the statistical 
relation is 

GM = H*(u,) (L ((,) H^(o,) (4-6) 

q 's 

where the (*) denotes the complex conjugate. All the other trans- 
formations implied by the matrix products in equation (4-4) are 
real, and they will eacfi have a power spectrum relation similar 
to equation (4-5) . The overall power spectrum relation between the 
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input vector process f and the output process vector y is then 

Gy = U H*M u"g^U M H (4-7) 

where the indication of the dependence on the frequency w has been 
dropped and the substitutions and = ri have been 

made. From an operational point of view the response spectrum 
matrix in equation (4-7) is obtained in the same manner as in equa- 
tion (4-6): that is, by premultiplying the input spectrum matrix 

by the conjugate transfer matrix and post-multiplying by the 
transposed transfer matrix hJ. 

The above procedure of first developing the vector signal 
block diagram and then drawing out the corresponding power rela- 
tions provides a powerful and convenient operational tool for 
more complicated processes. As a second example, consider the 
response of a continuous structure to a base input which is uniform 
spatially but random in time. The vector signal relation between 
a, the Fourier transform of the base input acceleration vector 
x(t), and the Fourier transformed absolute displacement response 
vector y has the block diagram representation shown in Figure 4-2. 
Taken as a whole, the transformations in Figure 4-2 constitute 




1 - T I .. 

+ — - I) M u' M^v 

OJ 


Fig. 4-2 Base Inp'Ut Transfer Matrix 
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(4-8) 


the physical-space transfer matrix 


H = U(H + ^ 

d L. 

( V ) 


1) M 


U^M 


To draw out the power relations for this process the same opera- 
tional approach used in the previous example is followed: 

= u(if + ^ 1) M iJm m u m \tt + I)LI^ 

(i) 


(4-9) 


Ul 


whore aqain matrix transposition lias been omitted for the synmetric 
-1 


matrices M, M , H. 

As a third example, the previous two examples are combined 
by assiminq that the force excitation and base excitation processes 
act simul taneously on the structure, and further, that they are 
statistically coupled. The vector signal block diagram for this 
case is shown in figure 4-J. 



tig. 4-J Combined Input Vector frocesses 

To draw ('Ut the power relations for the process shown in 
rigure 4-J tlie same operatiiMtal approach i- followed a^ with the 
previous examples. There is, howovor, a now Moik di.Hc.ur opera- 
tion of sunining s ta 1 1 s t ica 1 ly coupli'd vectors r.s’uices ,i 


n 




corresponding new power relation. This relation can be written In 
the following partitioned matrix form: 


-H, 




( 4 - 10 ) 


where 6^^ Is the cross power spectral density matrix which des- 
cribes the statistical coupling between the f and a Input pro- 
cesses. The compact relation of equation (4-10) Is more easily 
Interpreted In the expanded form; 

= "f\ "f - "f* ' "a* «af * "a* "a' > 

The first term on the right-hand side of equation (4-11) Is given by 
equation (4-7); the last term Is given by equation (4-9). Each of the 
middle two terms, because they represent cross power terms between 
two separate vector processes, are non-symmetr i c matrices. They 
are, however, conjugate transposes ot each other, so that their 
sum yields a matrix with the Hermitlan symmetry required for Gy. 

As a final example, consider the Input Identification problem 
associated with the system shown in Figure 4-3. It is assumed that 
the output process and only one of the two statistically coupled 
input processes are known. For example. Imagine that the signal 
vectors y(t) and f(t) are known and are available for statistical 
processing but that the spectral matrices and G^^ associated 
with the base input process are not known. If the cross process 
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matrix were known, It would In principle be possible to rearrange 
equation (4-11) and to reconstruct the unknown matrix G.. As an alter 
native, since y(t) and f(t) are available, it is possible to compute 
another cross power matrix G^^. The additional Information provided 
by G^^ makes possible the deduction of both G^ and G^^. We consider 
here only the reconstruction of 6_. The strategy for this approach 

Q 

is based on the following partitioned cross power relation for the 


system shown in Figure 4-3: 




The upper matrix equation of equation (4-12) is 


S ' - '■■fa 


(4-12) 


(4-13) 


Premultiply equation (4-13) by and substitute into equation (4-11) 


to obtain 


"a* =a ^ - <Nf* 'fy * 'yf 

This equation can be used to determine G. in terms of spectrum 

a 

matrices derivable from the input processes f(t) and the output 
process y(t). 

* T 

In a similar fashion an equation for can be 

obtained: 


G 


^H^ =G+H GH +(H G +G h') 
f ”f ”y "a a a ' a ay ya a ^ 


(4-15) 


This equation can be used to determine G^ in terms of spectrum 
matrices derivable from the input process x and the output process y. 



Response estimation by the use of equation (4-11) and Input 
Identification by the use of equations (4-14) and (4-15) appears 
deceptively simple. In fact, the amount of computation and 
measurement effort Implied by these equations can be quite large. 
The questions of computational strategy and efficiency are taken 
up in the next section. 

4.2.2 Computational Strategies and Efficiency: 

Program RAND is designed to efficiently compute power 
spectral density (PSD) and intearated power spectral density 
(IPSD) for two modes of operation: a forward (I.e., response) 
analysis and a backward (i.e., Input identification) analysis. 

A basic program strategy is that the major portion of the compu- 
tations are performed in a reduced modal space. The modal space 
is determined by examining the modal power participation 
matrices (described below) and selecting only those modes which 
actively generate significant response power. In general, the 
number of active modes m is less than the number of original 
physical degrees of freedom n, and the original system equations 
are reduced to a smaller and simpler uncoupled set. 

Further program efficiency is obtained by segregating the 
real and imaginary parts of complex matrices In such a way that 
complex matrix multiplication is replaced by a real dyadic 
operation called an overlay multiplication. Another program 
eff"'iency is provided which allows the user to spatially window 
a specific node in the system and to determine the power spectrum 
or integrated power spectrum for that node alone in all forward 
or backward modes of operation, thereby reducing computation as 
well as printout. 
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4.2.3 Method of Response Computation; 

Equation (4-11) can be recast Into a compact notation as 


G ® G G 1 • G ix 
r r/f r/a r/fa 


(4-16) 


where the "r“ subscript Implies the response variable which can 
be displacement (y), velocity (y), and acceleration (y). The 
symbol denotes "the response with respect to"; f implies force 
input, a the acceleration input, and fa, the cross power between 
the force and base inputs. Therefore, 6^^^ Is the partial res- 
ponse PSD due to a force input only; 6^^^^ is the partial response 
PSD for base acceleration only, and Is the partial response 

due to the cross PSD between the two random input processes. 6^ 
is the total response power and could be written as 


^r "" ^r/f, a, fa 


(4-17) 


As indicated earlier, the majority of the operations are performed 
in modal space. The transformation from modal space to the physical 
domain is performed by the operation 



(4-18) 


where the bar over the response PSD denotes a modal quantity. 
Thereiore, equation (4-16) can be reformed into a modal expression 
as 


^ ^r/f ^r/a " ^r/fa 


(4-19) 


The input power spectra are assumed to be separable in time and 
space; therefore, the partial power spectra become 


G = G(u)) D 


(4-20) 
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where 0 Is a matrix which describes the spatial correlation between 
the response degrees of freedom or the two random Input processes. 
The modal partial PSD In equation (4-19) can be written as 


5^/,. ' QfM Hf* R/ 

®r/fa • 


(4-21 ) 
(4-22) 
(4-23) 


where and are the coincident and quadrature components of 
the cross PSD of the force and base Inputs, * H Is the modal 
transfer function relating modal force to displacement, and 

= H + (1/u) ) I is the base acceleration modal transfer function. 

The 0 matrices are modal power participation matrices which 
are used in the response analysis to identify modes which actively 
generate significant power. The modal distribution matrices can 
be expressed as 


= M U' U M 

(4-24) 

-1 T 


= R U‘ M D M U R 

u a 

(4-25) 

-1 


Bfa “ U H U H 

(4-26) 


where the unbarred D's represent the spatial distribution matrices 
in the physical domain. Note that the modal power participation 
matrices need only be computed once and that they remain constant 
throughout the problem solution. Also, the complex operations In 
the above equations are performed using the special overlay 
multipl ications. 
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The complete output information for every degree of freedom 
is not always needed. In that case the response matrix can be 
spatially windowed at a selected node 1. This process will produce 
only the real co-power term since the quad-power term for a main 
diagonal element is always zero. A row matrix can be formed by 
partitioning the i^^ row from the modal matrix and designating 
it as pi. The scalar power at node 1 then becomes 

g’ • pT p, (4-27) 

Where the superscript i indicates that a windowing process has 
been imposed on the physical PSD matrix. 

The integrated PSD is obtained by an Euler integration pro- 
cess applied to matrices in modal space. This scheme was selected 
since it only requires the PSD matrices at step thereby avoiding 
the difficulties of saving previously evaluated matrices for inte- 
gration. Therefore, the integrated PSD is 

G(u). ) 

G{uj|^) = G(uj^_l) + — 2 — t^***!^ ^k-1^ (4-28) 

wliLM'e (') symbol denotes a nwdal integrated PSD and where 
- ^k^l ” new step, and ~ <i)|^ - is the 

old step. 

The rate at which the program steps through frequency solu- 
tion is detennined automatical ly by the program which uses 
the half power width of the nearest mode as its basic step in 


VJ 


frequency. In order to ensure convergence or to Improve accuracy* 
the program allows for subdividing steps and/or concentrating 
computations at selected modes of Importance. 

4.2.4 Method of Input Identification Computations: 

The general problem of identifying a spatially non-uniform 
input process is very complex both computationally and experi- 
mentally. A reasonable level of complexity can be obtained by 
assuming both force and base acceleration Input distributions 
are spatially uniform. Therefore, and are set to a one 
matrix and the input identification process is performed at a 
single nodal point on the structure which would correspond to 
an experimental accelerometer output. 

To estimate the base acceleration input PSD, equation (4-14) 
can be written as 


%/a ~ %/f,a,fa ^ %/f ~ ^fr/f,a,fa 


(4-29) 


wlicre Gj- is the partial power response due to the input base 
acceleration, to be identified, 

G is the total response of the system due to 

all effects. 

G : is the partial response due to the uniform force 
input rso, 

G.- ^ is the cross power between the input force 

^ PSD and the response of the system due all 
effects. 

As discussed previously, it is not reasonable to require an ex- 
periment to pix>duce all the data required to specify equation (4-29) 
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completely. Therefore, the windowing process must be applied to 
reduce equation (4-29) to 



r 


^ I O < 

V/f • V 


( 4 - 30 ) 


Note that the subscripts to the first and third terms have been 
dropped since they are the responses due to all effects. Equation 
(4-30) implies the following; 


g’ is the accelerometer response at nodal position i 

g’ is the input force PSD at nodal position i which 
' by assumption is the same at any node 

gI js the cross PSD information between the input 
force PSD and the output accelerometer at node i. 

This PSD contains a coincident and quadrature component. 


To detenitino an estimate of the input base acceleration PSD, 
the program applies a white noise PSD to the base in order to 
characterize the system power response, denoted as . Then 

the base input to node i can be written as 


g’ , = g’ (u) g’ , , 
r/a a^ ' r/a' 


(4-31) 


where G (w) is the identified base input acceleration PSD obtained 

ci 

from infoniiation at node i. Solving for g] using equations (4-30) 


and (4-31) produces 
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(4-32) 
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In a similar fashion, equation (4-15) can be windowed at node 


1 as 


^r/f * ^‘r ^’r/a ^ar 


i4-33) 
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where g!. Is the cross PSO between the base acceleration and 
accelerometer output as node 1. As In the previous case, a white 
noise force PSO Is applied to obtain the force power response 
^r/f • * identified force Input PSD obtained from Infor- 

mation at node 1 Is 

g{(.) . ^ IG^ ♦ gJ/. * g]^1 (4-34) 

Vf‘ 

4.3 Results of RAND Computations 

4.3.1 Combined Base and Acoustic Excitation of 
a Simple Oscillator: 

The combined base and acoustic excitation of a simple oscil- 
lator has been studied both analytically and by using Program 
RAND. Originally the exact analytical solution for this case was 
developed to rigorously check the program logic of RAND. However, 
since the simple oscillator is a prototype for each of 
the modes of a MOOF system, the analytical solution gives valuable 
insight into the convergence and resolution properties of the 
discrete frequency operations of RAND. 

For example, consider the acceleration response of a one DOF 
system of mass m, natural f'*equency and damping ratio c. For 
the combined force and base excitations characterized by the 
spectra CyCu), G (o), and the response Is 

G"{u)) = 



(1 G (uj) 


(4-35) 


where s ■ u/<d^. The term which contains the coincident canpo* 
nent of the cross power, appears with a negative sign. This reflects 
the fact that the separate responses to the force and to the base 
Inputs are In opposition to each other. In contrast, the term 
which contains the quadrature component of the cross power, 
appears with a positive sign and has a smaller magnitude than the 
term, at least for 2ce < 1. For the case where the Input spectra 
are constant with frequency the Integrated mean squared acceleration 
response Is 
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(4-36) 


The terms Involving v represent the overall output power contributed 
by the resonant portion of the transfer functions. Additional terms 
Involving & are non-resonant contributions to the output v4)1ch con- 
tinue to grow with frequency. Since the resonant contributions 

occur Just In the narrow band of frequency at w , a good estimate 

n 

of the Integrated response over a finite bandwidth which Includes 
the resonance can be obtained by dropping the limiting operation 
In equation (4-36). 

An Idea of the extreme range of the Input coupling effect can 

be obtained by examining the special case where ‘ 0 and 
2 

Gj » G^/m « C^^/m » Gq. This corresponds to the case where tl»e 
force and base inputs are fully coherent and the resonant contrl- 
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buttons of the two Input spectra are nearly equal. Substitution Into 
equation. (4-36) yields, for the finite upper frequency ratio 6, 

<y^>='6o«„B (4-37) 

Thus the resonant contributions cancel each other and only the 
non-resonant response remains. By changing the phase of the co> 
power term, = -6 q, the overall response becomes 

< = Vn 

where the power contribution due to resonance Is four times that 
caused by 6^ or alone. Figure 4-4 shows the results of a RAND 
computation for a related extreme case where f 0* The phase 
between the base and force input Is described by the phase angle 
~ tan 

In summary, the resonant power response ranges from near zero 
(coherent signals cancel), up through a factor of two (simple sum 
of uncorrelated powers) to a maximum of four (coherent signals add 
for a signal factor of two, a power factor of four). Note that 
this range will be reduced whenever the Inputs are not equal or 
are not fully coherent. It is clear that the effect of coherence 
and phase of the cross power can be very significant. Similar 
effects are apparent in MDOF system responses, to be discussed later. 

4.3.2 Acoustic Excitation of a Beam 

A simple ten degree-of-freedom beam was chosen to investigate 
the effects of the spatial correlation of the acoustic Input and 
the effects of the beam boundary conditions. Figure 4-5 shows the 
center node acceleration response PSD when the beam is clamped 
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Fig. 4-5 Effect of Pressure Correlation on Beam Response PSD 
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at each end and Is subjected to a base Input as well as either a 
spatially correlated pressure Input (solid line), or a spatially 
uncorrelated pressure Input (dashed line). All Inputs are cons- 
tant with frequency (white noise) and no correlation exists 
between the base or acoustic Inputs. It Is evident that the 
asymetrlc inodes (every other mode) do not participate In the 
response when the pressure Is spatially correlated. Participa- 
tion Is strong for these inodes, however, when the pressure Is 
spatially uncorrelated. The area under the uncorrelated curve 
Is, therefore, larger than that under the uncorrelated curve, 
resulting In a higher Integrated response for the uncorrelated 
case. 

Figure 4-6 shows the effect of boundary conditions on the Inte- 
grated acceleration power response. The solid line represents 
simple supports while the dashed line represents clamped supports. 
The discrete jumps In the curves occur every time that a partici- 
pating mode Is passed. Since the simply supported beam has the 
lowest natural frequency, its response starts to build first. 

The two curves then alternate in height for a time until the 
clamped, curve finally begins to maintain a higher Integrated 
value than the simply supported curve. These results Indicate 
that stiffen boundary conditions will give a somewhat higher 
Integrated response, but It Is felt that this will not be a major 
variable in response or input prediction. 

Figure 4-7 shows a family of backbone curves for Integrated 
response of the beam center node versus phase angle t|) = tan'^(Co/Qd) 
of the cross power between the base and spatially uniform acoustic 
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Inputs, for two support conditions. Omitted is the cos1ne*11ke 
variation with respect to each straight backbone line. (Each of 
actual curves In the family would have shapes similar to the 
curve shown In Figure 4-4.) Figure 4-7 demonstrates the substantial 
effect that statistical coupling between base and acoustic Inputs 
can have on the random response of structures. Similar effects 
will be demonstrated for a plate (Section 4.3.4). 

4.3.3 Single Point Force Excitation of a Plate 

During the development of RAND a number of solutions were 

run and comparisons made with existing analytical solutions In 

order to verify Program RAND. For example, RAND results were 

(3l 

compared with an analytical solution by Wlttlg- ' for the mean 
square velocity of a rectangular, simply supported plate subjected 
to a random point force having a white noise Input spectrum. The 
results of this comparison are shown In Figure 4-8. The RAND 
solution utilizes only sixteen degrees of freedom and therefore 
accounts for a very small portion of the inodes Included In the 
analytical solution. Still, however, convergence Is rapid and 
the two solutions are in good agreement. The question of conver- 
gence Is discussed further In the next section. 

4.3.4 Combined Base and Acoustic Excitation of a Plate 

Next we consider In more detail a plate response problem 

where the excitation Is more realistic than In the previous 
example. For this example an aluminum plate Is chosen which measures 
4 ft. by 3 ft., with a thickness of 0.2 in. This thickness Is 
chosen to Insure a fundamental frequency low enough to be affected 
by the low end of the Input spectra. The plate is assumed to be 
simply supported on all four edges. Kodal damping Is assigned a 
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value of c ■ .0& for the first mode and assumed to decrease for 
higher modes so that the product cu^ remains constant. Such 
a plate has s1xty>four natural frequencies below 1000 Hz. Only 
eighteen of these modes participate when excited by either a 
uniform pressure or a base Input. The Input spectra ranges In 
frequency from % to 1000 Hz as shown In Figure 4-9. 

At first the base and pressure Inputs are assumed to be 
statistically Independent. The Integrated displacement responses 
for these Independent Inputs, computed separately by RAND, are 
shown In Figure 4-10. Evident from this figure Is the fact that 
the displacement response converges rapidly for both Inputs, and 
Is almost entirely determined by the first few participating 
modes. One can conclude from this that any Instrumentation (such 
as a strain gage) that relies solely on displacement response for 
Input Identification would be completely Ineffective at mid and 
high frequencies. Because both Inputs are spatially uniform, 
the response Is symmetric and the largest response occurs at 
the center of the plate. The magnitude of the Integrated RMS 
response due to both Inputs Is 0.115 In., or approximately one- 
half of the plate thickness. 

The Integrated acceleration responses for Independent pressure 
and base Inputs are shown In Figure 4-11. They exhibit an entirely 
different character than the displacement responses In that con- 
vergence Is not evident. The Integrated response Increases a 
little as each participating mode Is encountered. Acceleration 
response data Is clearly preferable to displacement response data 
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Fig. 4-9 Combined Base and Acoustic Excitations for Plate 
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Statistically Independent Base and Acoustic Inputs 
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for the Identification of Inputs over a wide frequency range. The 
spatial distribution of acceleration over the plate (not shown 
here) Is again synmetrlCt and Indicates that the center of the 
plate has a higher response by a factor of almost two over any 
other plate location. This suggests that the center of the plate 
acts as a collector of energy and Is, therefore, the best site for 
the location of an accelerometer. 

The effect of statistical coupling between the base and 
pressure Inputs can be Investigated using the same plate model 
described above. The base and pressure PSD characteristics 
are chosen to be white noise, their respective power amplitudes 
sized to cause a roughly comparable response of the plate If 
applied separately. Fig. 4-12 shows the Integrated acceleration 
response at the center node of the plate for several variations 
In cross power phase angle v and for two values of the coherence 
ratio (all assuned constant with frequency). Note that the Inte- 
grated response can range from a low value of 4 g's RMS to a high 
of 30 g's RMS depending on the phase between the base and acoustic 
Inputs, a result consistent with the simple oscillator study 
discussed earlier In Section 4.3.1 For a phase angle of 90^ the 
response Is equal to 22 g's RMS. This value Is nearly equal to 
the response that would occur If the two Inputs were uncorrelated. 
This large range of possible responses Is directly attributable to 
the statistical coupling between the base and pressure Inputs, 
since these Inputs have remained unchanged In this plate response 
exanq)le. This emphasizes the need to assess the magnitude and 
phase characteristics cf any statistical coupling between Inputs. 
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4.4 Parametric Design Guidelines for Plates Subjected 
to Random Excitation 

Parametric relationships for plates subjected to random exci- 
tation are developed separately for a pressure Input and for a 
base Input. Consider first the forced response of a damped single 
degree of freedom system of mass N, resonance quality factor Q 
and natural frequency Proportionality relations for the mean 
square responses for displacement and acceleration are: 


6.Q 

mean square displacement - ^3 

M (i) 

n 

(4-39) 

mean square acceleration — ®f^n 

”17" 

(4-40) 


n 


where Is the Input force power spectral density taken to be a 
constant near The response of the single degree of freedom 
system can serve as a model for the resonant modes of a plate with 
the exception that the effect of a modal participation factor must 
be Included for each mode. Consider next that the effect of 
damping and modal participation factor are held constant, and only 
the plate mass and its natural mode frequencies are varied. The 
plate mass m is simply 

m = pabt (4-41 ) 

where p = mass density 

a,b - plate length and width 
t = plate thickness 
The fundamental frequency of the plate is 



2 D 

= “T 

n pt 
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(4-42) 


where 0 « plate bending stiffness 


. Et'’ p 


As the mass and fundamental frequency of the plate Is varied 
so will the mass and the frequency of each higher mode be varied 
accordingly. The total response to all the modes will then be a 
summation of modal responses, each similar to eqn (4-39) for dis- 
placement, or eqn (4-40) for acceleration. Parametric relations 
for the RMS response of the plate are obtained by substituting 
eqns (4-41) and (4-42) Into eqns (4-39) and (4-40), and then taking 
the square root of the resulting expression. These relations are 
shown In Table 4-1. The parameters can be divided Into three 
groups representing the effects of plate geometry, plate material, 
and plate thickness, (n refers to the material loss factor). 

These groups are represented by the symbols 0 and A for the displace- 
ment and acceleration response respectively, and with the subscripts 
6, M, and t to denote geometry, material and thickness. Figures 
4-13 through 4-17 present plots of the effects of geometry, material 
and thickness on plate response. 

Figures 4-13 and 4-14 show the dependence of the geometric para- 
meters Ag and Dg on the aspect ratio e of the plate for several 
values of the plate width b. Displacement parameter Dg decreases 
monotonically with increasing e. In contrast, acceleration para- 
meter Ag reaches a maximum at e = 1.33 and then slowly decays. 

The effect of material selection on displacement and acceleration 
response is shown in Figure 4-15. The displacement and acceleration 
parameters D|^| and A^| each shows a dependence on modulus E and 
density p which differs greatly from the modulus-to-density ratio 
commonly used in the design of precision aerospace instruments. 

As a result such inexpensive metals as steel and brass 
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Fig. 4-15 Effect of Material on RMS Plate Response 
(Pressure Input) 




compare favorably with the expensive metal beryllium as far as 
random response Is concerned. 

Figures 4-16 and 4-17 show the effect of thickness on response. 
The effect of thickness on displacement response Is very pro- 
nounced. In contrast, the effect of thickness on the accelera- 
tion response Is relatively mild. 

All of the parameter effects described above have been veri- 
fied by RAND. Typically the response of a particular plate can 
be determined for a given Input. The geometry, thickness and ma- 
erlal properties can then be changed arbitrarily and the new 
response predicted to within a few percent by using the parametric 
relationships. 

The development of corresponding parametric relationships 
for base excitation follow the same approach as for the pressure 
case. In this development there Is one major difference: It 
becomes evident that one can address either the relative or the 
absolute plate response. RAND studies show that the majority of 
the response of a plate to a base Input Is produced by motion of 
the plate relative to Its supports. Using D'Alembert's Principle, 
the forcing function for this relative response is proportional 
to the mass per unit area pt. This suggests that the parameters 
for base response can be obtained from the pressure response 
parameters simply by multiplication of the latter by the factor 
pt, as shown In the bottom half of Table 4-1. Note that only 
the geometry parameter Is unchanged. Changes In the material and 
thickness parameters will yield response predictions for the base 
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Input case which differ from the corresponding prediction for 
pressure input. As before, the accuracy of these base response 
parameters have been verified by Program RAND. 

The purpose of this parametric study was to develop the 
ability to predict changes in response due to changes in struc- 
tural configuration and as a result to be able to predict the 
response of a newly designed plate by comparing its physical 
characteristics to a plate with a known response. The above 
parametric relationships accomplish this goal. 

4.5 Conclusion to Part I 

Program RAND is a practical and useful tool for making res- 
ponse and input identification estimates for any multi-degree- 
of-freedoro structure subjected to botSi base and acoustic inputs. 
Displacement response computations, suitable for stress and sway 
space estimates, converge quickly and require only a low frequency 
representation of the structure. On the other hand, accelera- 
tion response computations, suitable for the estimation of g loads 
and for designing subassembly tests, do not converge quickly. 

In the latter case, the structural model must have good fidelity 
over a bandvndth which includes all frequencies of significant 
excitation. In the case of an acoustic input this bandwidth may 
extent to 10 KHz. 

Practical limitations to the RAND estimation procedure do 
not relate to the program itself, but rather to the program inputs. 
The limited fidelity of the finite element model used to charac- 
terize the structure has already been mentioned. A second input 
limitation, looking ahead to Part II of this study, results from 
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the difficulty of accurately measuring the random base or acoustic 
excitations. As far as model fidelity is concerned, perhaps a 
reasonable approach Is to use Program RAND for low frequency dis- 
placement studies, and for acceleration studies up to mid frequency 
range. Truly high frequency problems will probably best be studied 
using the technique of Statistical Energy Analysis. Input accuracy 
limitations are discussed in part II. 
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5.0 Part II: Exper1ment>Based Models for Structural Response 

mS- WdIR I3ehtTfk^!oh~ “ — 

5.1 Introduction to Part II 

In describing the experimental portion of this study a distinc- 
tion has been «iade between low resolution and high resolution ex- 
periments. This distinction reflects the historical fact that 
originally only a low resolution Fast Fourier Transform (FFT) ana- 
lyzer was available for use In this study. Subsequently a high 
resolution capability was added to the analyzer. Although the low 
resolution results are In some respects Inferior In quality to the 
high resolution results, there are advantages and disadvantages 
associated with both types of measurements. Consequently there 
Is much to be learned by contrasting the low and high resolution 
measurements, and they have been given equal emphasis In this 
report. 

5.2 Laboratory Test Description 

A simple schematic of the test set up for combined base and 
acoustic Inputs Is shown In Figure 5-1. The corresponding Instru- 
mentation layout Is shown In a more detailed fashion In Figure 
5-2. With this equipment it was possible to excite the plate 
by a base Input, by an acoustic Input, or by combined base and 
acoustic Inputs- Note that In the combined Input case both the 
base and acoustic excitation are caused by the same source, 
namely the MBT495 Random Noise Generator. This means, of course, 
that thc' base and acoustic Inputs are fully correlated statis- 
tically, and the nature of the cross power spectral density between 
these excitations is wholly determined by the frequency transfer 
characteristics of the base and acoustic channels which exist 
between the noise source and the aluminum plate. 





SUSPENDED ACOUSTIC CEILING 



Fig. 5-1 Test Set Up: Combined Base and Acoustic Input 
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In the experiments the procedure was to measure Input and res> 
ponse spectra under various conditions of excitation and then to 
compare these spectra with those predicted by Program RAND. In 
order for RAND to be executed It Is necessary to specify the dynamic 
modal parameters of the plate. The required parameters are the 
modal damping ratio vector the modal frequency vector f and the 
mode shape matrix U. Modal damping ratios can only be determined by 
experiment, and much effort was expended In this study attempting 
to obtain reliable damping estimates. On the other hand, frequency 
and mode shape estimates can either be made experimentally or by 
the use of a computer-aided finite element analysis. Probably the 
preferred approach is to determine these parameters wholly by 
experiment, since no idealized assumptions are required as with 
the finite element modelling process. And with the availability 
of computer-based modal analysis systems the experimental approach 
may also be easier. Unfortunately In this study no such modal 
analysis system was available. Mode shape estimates were therefore 
based on a finite element analysis using Program SAP IV. Damping 
and frequency data was obtained experimentally. 

The question of boundary conditions requires special attention. 
The analytical approach embodied in Program RAND utilizes the prin- 
ciple of linear superposition, and the ideal assumption is made that 
a base input, of infinite source impedance, acts independently of 
the acoustic input, which has zero source Impedance. The experi- 
mental reality is rather more complex. As will be detailed below 
in the treatment of input loading effects, both the shaker and the 
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loudspeaker possess a finite source Impedance. Further, these 
Inputs are dynamically coupled: for Instance, the base Input level, 
as measured by the plate fixture accelerometer will be caused to 
change merely by turning on the acoustic input. The RAND assumptions 
are still valid and the RAND predictive procedure is still appli- 
cable to real structures; It Is Just that particular care must be 
taken concerning boundary conditions In a test situation. Consider 
for example that we want to compare RAND predictions with experiment 
for the case of acoustic Input only. Since the analytical approach 
assumes infinite source Impedance for the base Input and, in this 
case, the base Input amplitude Is assumed to be zero, the proper 
experimental procedure Is clearly to mechanically block the shaker 
Input. Alternatively, a pure base input test, without any acoustic 
interaction, would require removing the loudspeaker, or, better yet, 
to remove any acoustic effects by operating in a vacuum. The diffi- 
culty with boundary conditions occurs only when It Is desired to 
test one Input at a time. Practically speaking, there will always 
be combined inputs In a field situation, and the RAND procedure 
will rigorously apply. 

Finally, we consider the question of data averaging. Whereas 
the analytically derived spectra produced by the RAND predictive 
procedure have in effect already been averaged mathematically, the 
experimental data produced by the base input and plate response 
accelerometers and by the microphone vary randomly, and must be 
sufficiently averaged to reduce the amount of statistical variation 
to acceptably small levels. A figure of merit for the adequacy of 
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this averaging Is the confidence level as measured by the equivalent 
degrees of freedom of the chi-square distribution model of the 
sampling process. Although the question of confidence level was 
addressed In these experiments It was not systematically monitored 
and Is not reported on here. A second figure of merit which Indi- 
cates the quality of data estimates Is the coherence ratio. This 
parameter Is related to the signal-to-nolse ratio of the data, and 
Is conveniently available on the FFT Analyzer. It was consistently 
monitored during all tests. A brief review of the coherence ratio 
follows. 

A common situation which arises In structural vibration measure- 
ment, which we shall refer to as output noise contamination. Is 
shown In Figure 5-3. 
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Figure 5-3 Output Noise Contamination 


In the Fourier domain, a linear dynamical structure, charac- 
terized by its transfer function H(«), receives an excitation X(u) 

and responds at a level W(w). An Independent noise source of strength 

2 

N(uj) Is present at the output. Define the coherence ratio y as 


2 

Y = 




(5-1) 
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where refers to the usual one sided power spectral density. 

Eqn (5-1) can be rewritten as: 

^ 

with 0<Y^<1 as 

Similarly, noise contamination Is possible at the input as 
shown in Figure 5-4. 


N(w) 



Figure 5-4 Input Noise Contamination 
Under the same assumptions, the coherence ratio is 

* r + 6""/6 

nn' XX 

2 

with 0<Y <1 as 

Both input and output noise contamination occurred during the experi- 
ments and will be discussed further below. 

5.3 Low Resolution Experiments 
5.3.1 Modal Damping Ratio Estimates 

Early damping measurements were made in the frequency domain 
using an FFT analyzer which possessed only a base band capability. 
Base band refers to the bandwidth of a spectrum which extends from 
zero frequency to a selectable upper frequency. The nominal 
frequency resolution of the analyzer is then the base bandwidth 
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divided by the number of spectral lines provided by the particular 
analyzer. configuration. Nominal resolution for the 400 spectral 
line analyzer (In base band operation) is given In Table 5-1 as 
a function of base bandwidth setting. 


Base Bandwidth 
Setti ng 

Nominal 

Frequency 

Resolution 

0 - 500 Hz 

1.25 Hz 

0 - 1000 Hz 

2.5 

Hz 

0 - 2000 Hz 

5 

Hz 


Table 5-1 Frequency Resolution for Baseband Operation 

As an example, low resolution estimates of the quality factor Q.j 

of the first plate mode (f^=240 Hz) varied from 25<Q^<100 depending 

on base band setting. The actual Q of the first mode, determined 

later by using an analyzer with a olgltal frequency expansion (Zoom) 

capability, was about 109. The actual half power bandwidth of the 

240 

first mode Is therefore about * 2.2 Hz. Clearly none of the 
above base bandwidth settings (Table 5-1) can adequately resolve 
the sharp resonance peak of the first plate mode. 

Time domain estimates of damping ratio were also attempted. 

The log decrement of the transient decay envelop of the plate were 
measured after passing the plate response signal through a band 
pass filter centered at the particular mode frequency of interest. 
Damping estimates were improved by this technique. For example, 
the first mode quality factor ranged from 76<Q<89. The method 
proved, however, to be tedious and inconvenient. It was also suspect 
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from a theoretical point of view: the plate dynamic characteristics 
were being windowed by the electronic filter which possessed Its 
own dynamic characteristics. 

The proper estimation of damping Is crucial to the success of 
any scheme of structural response or Input Identification estima- 
tion. And proper estimation of damping requires the use of Instru- 
mentation which has adequate resolution. It will be shown below 
In this report that there are certain advantages to making low 
resolution measurements of the Input and output spectra In structural 
dynamic studies. The requirement of high resolution measurements 
for damping, however, remains unchanged. 

5.3.2 Excitation and Response Spectra 

Typical examples of low resolution excitation and response 
spectra for the plate are shown in Figure 5-5, for the case of 
base Input, and In Figure 5-6, for the case of acoustic Input. 

In both cases a broad band Input excitation extends from nearly 
zero to 2000 Hertz. The five major peaks In the plate response 
spectra correspond to five participating normal modes predicted 
by a combined SAP IV and RAND analysis, as shown In Table 5-2. 


p 

SAP 


Experimental 

e 

Resonant 

RAND 

Resonant 

er 

Frequency 

Participation 

Frequency 


1 

252.2 Hz 

X 

250 Hz 

2 

391.1 



3 

622.8 



4 

626.7 

X 

605 

5 

755.5 



6 

954.6 



7 

981.2 



8 

1192 

X 

1105 

9 

1300 



10 

1324 



n 

1373 

X 

1310 

12 

1545 

X 

1445 

Table 
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Finite 

Element (SAP) 
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Fig. 5-6 Low Resolution Spectra: 


The SAP IV finite element model has yielded frequency estimates 
which are acceptably close to the experimental values, although 
the question of accuracy of the finite element model Is not an 
Important Issue, since we have already pointed out that the modal 
parameters (frequencies, damping ratios, and mode shapes) are 
best obtained experimentally using modal analysis equipment. Note 
that a spurious small response peak occurs at approximately 1200 
Hz. This was not predicted by the combined SAP/RANO analysis and 
was not Included as a transfer characteristic of the plate In any 
subsequent RAND response predictions or Input Identifications. 
Again, were modal analysis equipment available, this response 
detail would be routinely incorporated in the plate model repre- 
sentation of Program RANG. 

We turn now to the question of the quality of these spectra 
estimates. As stated above, the spectra have been sufficiently 
averaged timewise to the point of negligible statistical variation. 
However, a plot of the coherence ratio, shown In Figure S-7, Indi- 
cates that the signal-to-nolse ratio Is sharply reduced at a number 
of critical frequencies. In this Figure the base excited plate 
response Is presented as a semi log plot to allow study of regions 
of very large resonant response as well as regions of very low 
arti resonant response. Resonant and anti resonant responses are 
most easily described analytically In terms of the poles and zeros 
of the transfer function, which can In general be written In the 
form 


H(s) 


k(s-z^)(s-Z2) — 
(s-pj)(s-P2) — 



(5-4) 
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where the transfer function xeros z^, Z 2 *** are roots of the equation 
Z(s) = 0, and the transfer function poles p-ji P 2 *** are roots of the 
characteristic equation P(s) « 0. For very lightly damped structures 
these poles and zeros occur In conjugate complex pairs which are 
located very near to the Imaginary axis. The plate response spec- 
trum PSD Gp(w) is related to the base Input spectrum G^(u) by the 
equation 

fipW) - |H(»)|^ <5-5) 

where H(Ju) Is the frequency response function obtained from the 
transfer function H(s) simply by setting s>jw. A resonance occurs 
when |H(ju)| Is maximum. This In turn occurs when the value s«jw 
Is very close to (effectively Is equal to the Imaginary part of) a 
transfer function pole p^, P2~“~ Similarly an antiresonance 
occurs when |H(Jui)| Is a minimum. This In turn occurs when the 
value S-jui passes very close to (again, effectively equals the 
Imaginary part of) a transfer function zero z^, Z 2 — etc, A 
resonance Is a global characteristic of the structure and exists 
independently of the spacewlse nature of the Input. An anti resonance, 
on the other hand. Is a local characteristic of the structure; that 
Is, it applies to a particular location only. Further, the location 
of the anti resonance, both spatially and In frequency, depends on 
the spatial distribution of the Input excitation. 

Consider now the two antiresonant responses which occur approx- 
imately at 825 Hz and at 1350 Hz as shown In Figure 5-7. Since 
the damping of the structure Is very small, the antiresonant res- 
ponse Is also very small (zero damping would result In truly a zero 
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response). Clearly such a small signal will be corrupted by ever 
present backgri...id noise. Hence a sharply reduced value of the 
coherence ratio. Since this noise contamination occurs at the 
output (response) of the transfer function H(w) we have an example 
of output noise contamination as shown In Figure 5-3. 

' Significant reduction In the coherence ratio also occurs at 

the resonant frequencies. However, the mechanism for this - an 
Input noise contamination caused by Input loading effects - Is 
obscured in this data because of Inadequate resolution of the FFT. 
The details of the input loading phenomenon will be taken up later 
during the discussion of the high resolution experiments. 

5.3.3 Response Prediction Using Program RAND 
At this point an attempt was made to compare the pUte response 
level obtained experimentally with the same response level predicted 
by Program RAND. The RAND model of the plate used damping ratios 
determined by free vibration decay measurements, natural frequencies 
obtained by FFT measurement, and mode shapes obtained by finite 
element analysis. Four separate experiments were conducted using 
base input excitation only. The base Input and plate response 
spectra were measured using the low resolution FFT analyzer. Be- 
cause the FFT analyzer could not resolve the sharp response peaks 
adequately, there was little hope that experimentally measured 

^ response spectra would agree with the corresponding spectra pre- 

► 

^ dieted by RAND. However, FFT measurements of the rather slowly 

; varying broad band input spectra could be expected to be fairly 

good. And the integrated area of the response peaks, that Is, the 



overall power of each resonance* would also be fairly accurate 
because the FFT algorithm accurately computes the correct power of 
the time windowed signal* Independent of resolution adequacy. The 
results of these four runs are given In Table 5*3. The agreement 
between experiment and Program RANO Is rather satisfactory* parti* 
cularly when It is remembered that the excitation data Input to 
RANO suffered from noise contamination and t^e Input spectrum was 
poorly resolved. Reasons for the relative success of the predic* 
tion procedure In the presence of these error producing effects 
will be discussed later after the high resolution results have 
been presented. 


Pun 

No. 

Total Overall Response 
Experimental Program RANO 

Error 

1 

40.4 g^ 

44.8 g^ 

+.45 dB 

2 

41.1 

35.2 

-.67 

3 

54.3 

40 9 

-1.2 

4 

48.1 

37.3 

-1.1 


Table 5.3 Plate Response Estimates: Base Input, Low Resolution Spectra 

5.3.4 Input Identification Using Program RANO 

The Input identification mode of Program RANO takes as input 
data the response spectrum of a structure, accounts for the transfer 
characteristics of the structure, and then computes the spectrum 
of the excitation applied to the structure. The use of a low 
resolution FFT analyzer for input Identification mode measurements 
will not be satisfactory since the analyzer must resolve thj details 
of a response spectrum which consists of many sharp resonant peaks. 
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Nevertheless, It Is Instructive to look at an early example of an 
input spectrum Identification provided by the RAND procedure using 
data obtained by the low resolution FFT analyzer, shown in Figure 
5-8. A comparison viith the actual input spectrum (not shovm) 
shows that the input spectrum estimate is of poor quality as 
anticipated, particularly near the resonant frequencies at 250 Hz, 
605 Hz, 1105 Hz, 1310 Hz and 1445 Hz. But the most important aspect 
of this spectrum estimate is that it is dominated by three peaks 
which are very obviously in error, since the actual spectrum v«s 
relatively flat with frequency. The presence of these peaks has 
nothing to do with the low resolution of the FFT analyzer. The 
first two peaxs, at 825 Hz and at 1350 Hz, are caused by the poor 
signal -to-noise ratio at the two a.nti resonances described above 
in Section 5.3.2. At an antiresonant frequency the transfer 
characteristic of the structure decreases nearly to zero. The 
structural response, as measured by an accelerometer and estimated 
by an FFT analyzer, does not decrease nearly to zero ( as it should 
theoretically, and as Program RAND expects) because of the presence 
of noise. The only explanation for this non-vanishing response, 
as far as Program RAND is concerned, is that the input must have 
a very sharp peak. Thus erroneous peaks will always appear in 
the input identification estimate at each anti resonant frequency 
of the structure. From a practical point of view the anti resonance- 
induced error can be avoided by simply instructing Program RAND to 
set the input identification estimate to zero for a small band of 
frequency about each anti resonance, a process called blanking. 
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The implication of this phenomenon is that the transfer charac- 
teristic. of the structure which is stored in RAND must explicitly 
include antiresonant frequencies as well as resonant frequencies. 

It is also clear that the input identification procedure will 
break down in the region of an anti resonance, and that a portion 
of the input spectra will be lost because of an inadequate signal- 
to- noise ratio. 

The third large source of error in Figure 5-8 extends over a 
band of frequency 1800<f<2000 Hz. This is caused by the fact that 
the transfer characteristic stored in Program RAND includes only 
five modes and is therefore only representative of the structure 
out to a frequency of, say, 1600 Hz. The error can be removed by 
simply including an additional (sixth) mode, thereby extending 
the fidelity of the transfer characteristic to 2000 Hz. 

The problems described in this section were discovered during 
early low resolution input identification runs, and they called 
attention to the importance of assuring adequate bandwidth for 
the structural transfer characteristic and to the value of anti- 
resonant frequency blanking. These refinements were routinely 
incorporated in all subsequent low and high resolution studies. 

5.4 High Resolution Experiments 

5.4.1 Improved Estimates of Damping Ratio and Frequency 

The addition of high resolution capability to the FFT analyzer 
quickly led to damping ratio and resonant frequency estimates which 
were reproducible and also independent of any further reduction of 
the resolution bandwidth setting of the analyzer. With this capa- 
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blllty for precision It soon became apparent that the damping ratios 
and resonant frequencies were affected by subtleties in the boundary 
conditions of a particular test: for Instance, whether or not the 
shaker fixture was blocked, what type of excitation was used, and 
whether or not the loudspeaker was In position over the plate. The 
variation in frequency with changes In boundary condition were 
typically less than 1%. Unfortunately, due to Input loading effects 
such small variations In frequency are still significant as far as 
structural response and Input identification estimates were con- 
cerned. Discussion of the effect of boundary conditions will be 
postponed until after an examination of Input loading models has 
been made. 

Representative results for high resolution frequency and 
damping ratio are given In Table 5-4. The tabulated values are 
based on the base Input transfer function obtained by dividing the 
plate response (at center node) by the base Input (at plate fixture) 
Theoretically (and experimentally) these results are the same as 
for the boundary condition case In which the fixture motion Is 
blocked. 


Mode 

Number 

Resonant 

Frequency 

Damping 

Ratio 

Q 

1 

239.2 Hz 

.0046 

109 

2 

575.5 Hz 

.0022 

227 

3 

1050.3 Hz 

.0024 

208 

4 

1260.5 Hz 

.0015 

333 


Table 5-4 High Resolution Mode Frequencies and Damping Ratios 
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These damping ratios are significantly lower (and the corresponding 
Q values higher) than values obtained previously by low resolution 
FFT analysis or by the time domain transient decay technique. 

Once more we emphasize that adequate frequency resolution is essen- 
tial to the attainment of realistic damping estimates. 

5.4.2 Input Loading Effects: Base Excitation 
5.4. 2.1 First ilode Input and Response Characteristics 
High resolution dynamic characteristics of the first mode 
obtained from a random base excitation test are shown in Figure 
5-9. Several features are of Interest. First, the base input 
(curve a) shows a sharp decrease and subsequent increase In level 
as the plate response peak (curve b) Is passed. The Input minimum 
is substantial enough to cause the coherence ratio (curve c) to 
decrease, showing that the Input minimum Is contaminated by noise. 
Curve d shows the phase characteristic of the transfer function 
H(tu): specifically, the phase of the plate response lags behind 
the base input phase by 90° at a frequency of 238 Hz. Finally, 
the magnitude of the transfer function (curve e) reaches a maximum 
at the same frequency. Note that this frequency differs from the 
frequency at which the plate response magnitude is a maximum. 

This pattern, or a similar pattern, rather generally occurred at 
each structural resonance under both base excitation or acoustic 
excitation. The pattern consists of a decrease in the input level 
attended by a decrease in the coherence ratio, also an increase In 
the input level, and a response maximum which occurs at a frequency 
which differs from the frequency of maximum transfer function 
magnitude. 
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This Input loading effect Is not a new phenomenon. Much 
attention was given to Input loading In the early days of vibration 

4 

testing by the designers of electrodynamic shaker systems . In 
fact, since the loading phenomenon was considered an undesirable 
feature of vibration testing, electronic compensators called peak- 
notch filters were developed to render the vibration Input spec- 
trum flat with frequency. More extensive need for input spectrum 
adjustment led to the development of manual and finally automatic 
compensation devices known as equalizers. 

Our Interest here in the Input loading phenomenon stems from 
the expectation that Input loading should be a relatively common 
occurrence In field studies of real structures. The value of ana- 
lyzing the details of the Input loading phenomenon in the rela- 
tively simple laboratory set up of this study Is the Insight that 
It provides for the more complicated field system. 

5. 4. 2. 2 Model Development for Input Loading Analysis 
As stated earlier, the Input loading Is caused by the fact 
that the base and the acoustic Inputs possess finite source Impe- 
dance. Consequently, the base Input model roust Include a real 
displacement source of finite Impedance rather than an Ideal dis- 
placement source with Infinite Impedance. Similarly, the acoustic 
Input must be modelled as a real force source of finite Impedance 
rather than an Ideal force source of zero Impedance. A simple 
model which adequately explains the Input loading effects experi- 
mentally observed at the first plate resonance Is shown In Fig. 
5-10. Parameters used In Fig. 5-10 are defined as follows: 
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Fig. 5-10 Model for Input Loading Analysis 


2 



} 

\ 

I 

f magnetic force applied to speaker cone (zero 

° Input Impedance) 

m, mass of rigid speaker cone and exciter coll 

o 

k^,b^ stiffness, damping coefficients for cone suspension 
impedance of acoustic transmission line 

^ m^ plate mass associated with first plate mode 

k ,b stiffness, damping coefficients associated with 

P P first plate mode 

f magnetic force applied to shaker armature (zero 

input impedance) 

m^ rigid body mass of shaker armature and plate 

f 1 xture 

k ,b stiffness, damping coefficients of shaker armature 

^ suspension 

The analysis of the input loading phenomenon Is most easily 
presented as two special cases: base excitation of the plate. Ig- 
noring the acoustic elements, and the acoustic excitation of the 
plate with the shaker input set to zero (blocked). These cases are 
considered below in sections 5.4.2. 3 and 5. 4. 3. 2. 

5. 4. 2. 3 Base Input Loading Analysis 

An example of high resolution excitation and response spectra 
for a base input run was given in Section 5. 4. 2.1, as well as a 
broad description of the main features of the input loading effect. 
We proceed now to study the base input case in more detail. 
i Consider the case of base Input only. The acoustic elements 

' are assumed to be absent. The model of Fig. 5-10 simplifies to a 

two degree of freedom representation of the plate and shaker system. 
The plate system, by itself, has a natural frequency u and damping 

I * * 
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ratio Cp. The shaker system, by Itself, has a natural frequency 
and damping ratio If the resonant frequencies of the 
shaker anti of the plate are widely separated (wp»Wj), the fourth 
degree characteristic function which determines the coupled 
resonant frequencies of the combined system can be factored Into 
two quadratic terms as follows: 

characteristic function* (-C 4 )(-« (5-6) 


that Is, the characteristic frequencies of the combined systems can 
be described by a lower resonant frequency and a higher 

resonant frequency The damping ratios of the separate 

systems, the shaker damping ratio and the plate damping ratio Cp» 
are also affected when the shaker and plate are coupled. A first 
approximation to the relation between the separate and the coupled 
systems are: 



0 “c 

(1 . _E . 4r r -1) 
“s'* m^ “p 




c 

h 






^^s'^p «p^ 



(5-7) 

(5-8) 

(5-9) 

(5-10) 


The transfer functions for the base Input and plate response accele- 
rations can then be written as: 
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(“« + j2c <D (1) + u ) 

— 2 ' 2 “ — !T 

(-W +j2Cjj^Uj^u+(o^ )(-u w+«^ ) 

-u)^ fj/m^ 

— 2 — 2 2 — 

(-0J +j2c^u)j^u+u)j^ )(-<«» +j2C|jU^ui+u^ ) 


(5-11) 


(5-12) 


If now we evaluate these transfer functions near the plate resonant 
frequency WpW the following block diagram is obtained: 


- + j2fpWj, + 

Vs 

+ Wp2 


- 


- + iZfpWpCJ + Wp^ 



Figure 5-11 Base Excitation Block Diagram 

Several aspects of this diagram are of particular importance. First, 
the transfer function on the left consists of a conjugate pair of 
zeros and a conjugate pair of poles. On the complex plane, these 
roots are very close to each other and also very close to the ima- 
ginary axis, since and Wp differ only by very small quantities 
and the damping ratios and Cp are «1. The corresponding fre- 
quency response function has precisely the same anti resonant and 
resonant shape as base input curve (a) shown in Fig. 5-9. A second 
important aspect of the diagram is that the right hand transfer 
function precisely corresponds to the transfer function H(u) shown 
as curves (d) and (e) on Fig. 5-9. The final and perhaps most 
important aspect of the diagram is the fact that if these two 
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transfer functions are cascaded to form a single transfer function 
between the shaker input force and the plate response acceleration, 
there will be a pole*zero cancellation associated with the plate 
frequency tOp as shown In Fig. 5-12. This leads to the simple, 
and perhaps surprising result: since Wp « the transfer func- 
tion shown in Fig. 5-12 Is nearly that of a simple oscillator 
excited by a constant Input acceleration f^/m . 

i 5 


fj/m 


$ 



i2fpWpW ♦ Wp^ 






Figure 5-12 Shaker Force to Plate Acceleration Transfer Function 

This particular input acceleration level is the value that would 
occur if the plate were to cause no loading effect at all. In 
essence, then, the net effect of the input loading effect on the 
plate resonance is to shift it slightly upwards in frequency from 
ujj to 0 ^. The magnitude of the plate resonance is virtually un- 
changed from the value it would have in response to an unloaded input 
fj/m^ which is constant with frequency. So the effect of the input 
loading phenomenon on the plate response is rather minimal. How- 
ever, the effect of the input loading phen-^menon on the input to 
the plate remains problematic. As we have seen, the input loading 
effect suppresses the input sufficiently to allow the signal to be 
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i 

I 

I contaminated by noise. It also complicates the problem of Input 

i' identification because It Is necessary to accurately resolve all 

I 

the details of an Input which varies very rapidly with frequency. 

and to match the peak and notch frequencies of the Input curve 

very accurately with the plate response peak frequency. 

5. 4. 2. 4 Response and Input .’dentificatlon Estimates for 
' Base Excitation 

The use of more accurate mode damping ratios and frequencies 

led to modest improvements in response estimates and substantial 

improvement in input Identification estimates. The results of 

plate response experiments with base excitation are shown In 

Tables 5-5 and 5-6. 

Table 5-5 High Resolution Base Input Test; Response PSD at Resonance 




at Resonance 

Error 

Experimental 


1 

3.6 g^/Hz 

5.13 g^/Hz 

+1.5 dB 

2 

9.18 

7.14 

-1.1 

3 

4.69 

3.43 

-1.4 

4 

4.19 

3.02 

-1.4 


Table 5-6 High Resolution Base Input Test: Modal Response Power 


Mode 

No. 


Error 



1 

21.1 g^ 

26.3 g^ 

+ .96 dB 

2 

46.9 

38.7 

- .83 

3 

39.9 

29.1 

-1.4 

4 

39.3 

21.5 

-2.6 

Overall 

155.4 

127.9 

- .85 














1 


Table 5*5 shows the peak response PSO obtained by experiment and 
by RANG. These results represent an Improv^nt over the results 
attainable by low resolution experiments (low resolution measure- 
ments of the sharply peaked power response spectra were typically 
In error by 12 to 16 dB). On the other hand, discrepancies between 
the theoretical and experimental modal power estimates (Table 5-6) 
are approximately the same as for the low resolution estimates 
(Table 5-3). As mentioned before, this Is because the FFT algo- 
rithm computes the correct power of a signal Independent of the 
degree of frequency resolution. 

Substantial improvement In Input Identification accuracy Is 
obtained by using high resolution FFT measurements. This Improve- 
ment In accuracy Is particularly evident at the plate resonant 
frequencies where large errors had previously existed due to 
Inadequate resolution of input loading effects. In Figure 5-13, 
the RANG estimation of the base Input In the region near each of 
the four plate resonances Is compared with the measured Input data. 
The plot Is semi logarltimilc. By comparing both the measured and 
identified base Input power spectra to an arbitrary reference 
value, the ordinate can be scaled In dB as shown. The difference 
between the two curves can then be interpreted directly as an 
error expressed In dB. The rather coarse sampling of these two 
durves reflects the fact that the data link between the FFT analyzer 
and Program RANG was not automated and input/output data taking was 
done by hand. The accuracy of the input identification procedure, 
however. Is satisfactory. An automated data link would improve the 


73 




dB 


3 I 
2 ^ 

1 f 

Q ^ 



Fig. 5-13 High Resolution Identification of Base 
Input to Plate (Four Modes) 
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accuracy and resolution considerably. The measured overall power 

2 

for the frequency band 0<f<2000 Hz was .60 g . The corresponding 

2 

input identified by Program RAND was .524 g , which represents an 
overall error of -.59 dB. 

5.4.3 Input Loading Effects: Acoustic Excitation Experiments 
5.4. 3.1 First Mode Input and Response Characteristics 
An example of low resolution excitation and response data for 
an acoustic excitation test has already been given in Fig. 5-6. 

Two examples of high resolution data which show the effect of acous- 
tic input loading on the first plate resonance are given in Figures 
5-14 and 5-15. Broadly speaking, the acoustic loading effect is 
similar in pattern to the base input loading effect, although a 
closer study will reveal differences in detail. The actual shape 
of the loaded acoustic input, for instance, is different from the 
shape of the loaded base input curve. This is because the dynamics 
of the acoustic excitation differs from the base input dynamics. 
Dynamic models for the acoustic excitation of the plate are dis- 
cussed in the next section and in the Appendix. An interesting 
contrast also exists between Figures 5-14 and 5-15. The first 
figure, which refers to the test condition of a blocked fixture 
(refer to Figure 5-10; yjj=0), exhibits a relatively mild acoustic 
loading eff'^cf. The second figure, which refers to the test con- 
dition of an unblocked fixture, exhibits a relatively strong 
acoustic loading effect. (The relative strength of the loading 
effect is measured by the depth of the anti resonant region of the 
pressure curve and by the corresponding amount of input noise 


80 




243 Hr 


f 



Fig. 5“iti High Rasolution Spectra; Plate Response to 

Acoustic Input (First Mode, Fixture Not Blocked) 




contamination, as evidenced by a decrease in coherence). The ex- 
planation for this difference is that the method of fixture 
restraint used in the blocked fixture test caused the damping 
level to be larger than that for the unblocked fixture test. 

This emphasizes the fact that the input loading phenomenon becomes 
more evident as the structural damping level decreases. 

5.4. 3.2 Acoustic Input Loading Analysis 

Although the basic physical explanation for the acoustic input 
loading and the base input loading cases is similar, the acoustic 
and mechanical models required for the acoustic loading analysis 
is too complex to include in this section. These models are pre- 
sented in some detail in Appendix A. Two acoustic excitation 
cases are of interest: excitation from a distant sound source 

which is not loaded by the plate generated sound waves, and exci- 
tation from a nearby sound source which is affected by the presence 
of the plate. The later case represents the experimental set up 
used in this study. 

The results of the distant sound source analysis are presented 
in block diagram form. Fig, 5-16 (a). An incident pressure wave p^ 
emanates from the distant sound source. It is partially reflected 
by the moving plate. The resulting total pressure at the plate 
surface, p^, is the sum of the incident and reflected wave pressures. 


2 (- +j2fpajpU) + Wp^) 


CM 

3 

< 

-cj2+j2fp'ujpa) +Wp^) 




<a) 


2A 


co‘ 


-w^+i2fp-cOpW+Wp2 


(b) 


Figure 5-16 Acoustic Input Loading; Distant Source 


The transfer function on the left again consists of a conjugate 
pole-zero pair. Although the undamped resonant frequency of the 
pole and the zero pairs are the same, the pole pair is more 
heavily damped. Specifically, the pole damping ratio 


bp -t- pcA 


(5-13) 


2 V k m^ 
P P 


includes both the material damping coefficient bp of the plate and 
the acoustic radiation damping effect pcA. The zero pair damping, 
on the other hand, includes only the plate material damping coef- 
ficient bp. The frequency response characteristic of the left 
hand transfer function is a gain of 2 (perfect reflection) at all 
frequencies except at resonance, when it drops sharply down to 
t'.a value 2b /(b +pcA). Thus an input loading effect exists 

r r 
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even though the distant acoustic source is not affected by the 
reflected wave from the plate. Such a frequency response charac- 
teristic does not cause a shift in frequency of the plate accelera- 
tion response peak. The transfer function on the right is again 
simply the transfer function for the plate under acoustic exci- 
tation. When cascaded, as in Fig. 5-16 (b), the plate acceleration 
response exhibits a single resonance, with the same amplitude it 
would have if excited by the Incident pressure p^ without loading 
effect. 

Next we consider the case of a nearby sound source (See 
Appendix for details). In Fig. 5-17 (a), a magnetic force f. acts 

a 

upon the loudspeaker cone, causing a pressure p^ to act upon the 


f ^ 

-w^mp+jcobp+kp 

Po 


V 


A 


-w^mp + jwbp + kp 



(a) 



(b) 


Figure 5-17 Acoustic Input Loading: Near Source 

plate. The characteristic function A of the left hand transfer 
function is a complicated relation between the parameters which 
characterize the loudspeaker, the plate, and the connecting 
acoustic transmission line; 
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(5-14) 


L = f-(pcAw)^ + (-w^m +jcub +k. )(-u)^m +ju)b +k_)1 sin «L/c 

L uaa PPPJ 

+ pcAu) cos uL/c 

The resonant frequencies defined by the characteristic equation 
have not been evaluated explicitly, since a detailed evaluation of 
the many parameters would be required which is beyond the scope of 
this study. Note, however, that again the left hand transfer 
function which defines the pressure input to the plate consists of 
a conjugate pole-zero pair, and that again a pole-zero cancellation 
occurs when the two transfer functions are cascaded as shown in 
Fig. 5-17 (b). Our conclusion is again the same: although there 
is a complicated input loading interaction between the plate and 
the acoustically coupled loudspeaker system, the plate response 
is nearly the same as it would be if there were no input loading 
effect at all. The actual resonant frequency of the acoustically 
coupled plate differs slightly from the plate frequency which is 
excited by a distant pressure source. And again, the estimation 
of details of the acoustic input p^ remain problematic for reasons 
already stated above: the difficulties of estimating and measuring 

of high resolution narrow band spectra which vary rapidly with 
freqi' ncy. 

5. 4. 3. 3 Response Estimates for Acoustic Excitation 
The acoustic input and the acceleration response of the plate 
were measured and compared with the plate response predictions 
of Program RAND. Table 5-7 shows the resonant peak PSD estimate 
comparisons for the first four plate resonances. Table 5-8 shows 
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Table 5-7 High Resolution Acoustic Input Test: 
Response PSD at Resonance 


Response PSD 

at Resonance 

Experimental 

Program RAND 

.00748 g^/Hz 

.00616 g^/Hz 

.0283 

.0156 

.0300 

.0138 

,0240 

.0432 


Error 


-.087 

- 2.6 

-3.4 

+ 2.6 


Table 5-8 High Resolution Acou stic Input Test 
Modal Response Power 


Response 

Power 

Experimental 

Program RAND 

.025 g^ 

.024 g^ 

.084 

.085 

.234 

.175 

.123 

.503 


Error 













the integrated power comparison of each of the four resonances. 
These results are comparable to the plate response estimates ob- 
tained from base excitation tests. In both the base input and 
the acoustic input tests very close agreement between theory and 
experiment is possible. Yet it is also possible for substantial 
discrepancy to occur (up to 6 d6 in one acoustic input case), 
even in spite of the fact that measurements were made with a high 
resolution analyzer. These large errors result mostly from imper- 
fect characterizations of the sharply varying loaded inputs. 

In an attempt to avoid the difficulties of measuring the 
sharply varying acoustic input spectra, an alternative method of 
characterizing the input was introduced: the taking of a simple 

average of the input PSD over a frequency bandwidth centered at 
the plate resonance. This approach requires less effort than the 
point-by-point characterization of the input, since the average 
is simply the ratio of the integrated power (obtained from the FFT 
analyzer by either a high resolution or a low resolution measure- 
ment) divided by the bandwidth. The approach also yields better 
results, as shown in Tables 5-9 and 5-10. A theoretical basis 
for the relative success of the bandwidth average approach has 
already been developed in the input loading studies in Sections 
5. 4. 2. 3 and 5. 4. 3. 2. In both these analyses it was found that 
a sharply varying loaded input causes only a slight shift in the 
plate response resonance and does not significantly affect the 
response amplitude. It is reasonable that the loaded input can 
just as well be replaced by an equivalent input which is constant 
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Table 5-10 High Resolution Acoustic Input Test: Modal Power 
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with frequency. The easiest equivalent Is a simple bandwidth 
averaged Input. Further discussion of this simplified approach 
Is presented below In Section 5.6. 

5.5 Combined Base and Acoustic Excitation Experiments 

The base excitation and the acoustic excitation were simul- 
taneously applied to the plate at amplitudes which would each 
separately cause comparable response levels. Low resolution 
measurements were made of the pressure input PSD, the base accelera- 
tion PSD, the coincident and quadratu;« ci<nponents of the cross 
PSD for the pressure and base inputs, and the integrated response 
power of the plate. These spectra are shown in Figures 5-18 through 
5-20. No corresponding response estimates have been made using 
Program RAND. 

Two aspects of these spectra are of interest. First, each 
of these inputs, applied separately, would cause comparable res- 
ponses of the second, third and fourth modes. When these Inputs 
are combined, however, the plate response shows a relatively en- 
hanced second mode and a relatively attenuated third mode. This 
effect Is caused by the phase relationship between the two Inputs. 
Analysis of the response of a simple harmonic oscillator to 
statistically coupled inputs (Section 4.3.1) showed that the two 
inputs can either work together or work against each other, de- 
pending mostly on the sign and magnitude of the coincident com- 
ponent of the input cross power spectral density. A large negative 
coincident component will cause reinforcement of the inputs, whereas 
a large positive coincident component will tend to cause input 
cancellation. On this basis, the first mode should be suppressed, 
the second mode enhanced, the third mode unaffected, and the 
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Cross Spectra for Combined Force and Acoustic Inputs 







^SD for Combined Base 




I fourth mode strongly enhanced » each relative to the superposition 

of the response to the acoustic and base Inputs applied separately. 
Roughly speaking, these effects appear consistent with the experl- 
i mentally measured responses shown In Figure 5-20. 

I A second feature of the combined excitation spectra Is the 

‘ obvious presence of Input loading effects at the resonant fre- 

quencies of the plate. This Is particularly noticeable on the 
coincident and quadrature components of the Input cross power 
I (Figure 5-19) even though the low resolution FFT analysis would 

generally tend to obscure such effects. 

It has already been established that it is difficult to 
obtain good response estimates when they are based on measured 
! Input spectra which experience Input loading effects. The case 

t ■ 

will be no different as far as the input loaded cross power 
spectra are concerned. Although a bandwidth averaging technique 
offers some improvement for the case of individual input spectra, 

I . it remains to be demonstrated whether or not a similar procedure 

will be successful when cross input spectra are concerned, 
j 5.6 Simplified Treatment of Input Loading Effects 

I 

We have already seen the effectiveness of using a simple 
i bandwidth averaging technique when dealing with input loaded 

, spectra. This technique simplifies the measurement procedure 

; and leads to a better structural response estimate. Some ques- 

I 

tions remain, however. Analysis of the input loading phenomenon 

i 

i suggests that a constant, equivalent Input spectrum can replace 

I . ^ 

I 

the sharply varying input loaded spectrum. The analysis does 


"1 
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not suggest, however, that the bandwidth average is the correct 
estimate of the constant, equivalent Input. As damping decreases, 
for Instance, the bandwidth average Increases, whereas the theo- 
retical value of the constant, equivalent Input remains unchanged. 
There Is also the practical question of determining a suitable 
bandwidth for the average. A larger bandwidth lessens whatever 
error may apply because the technique Itself Is not theoretically 
correct. On the other hand, a large bandwidth will tend to obscure 
Input spectrum variations Which have nothing to do with the Input 
loading phenomenon. The answers to these questions will require 
further study. In the mean time, the bandwidth averaging technique 
offers a distinct Improvement over precise poInt-by-poInt measure- 
ment of the Input loaded spectrum. 

What procedure should be followed when Input Identification 
estimates are to be made? Not only Is such an Identification of 
an Input loaded spectrum difficult, but it would seem that It Is 
also undesirable. It Is of much more practical Interest to Iden- 
tify the constant, equivalent Input spectrum which will yield the 
same measured resonant response. In this case a theoretical guide- 
line Is available. For an Isolated resonant mode excited by a 
white noise the acceleration response power Is given by 


mode response power 



Go Q 

on 


(5-15) 


where J’ Is the modal participation factor. In the input Identifi- 
cation situation, the white noise Input G^, equivalent to the Input 
loaded spectrum near the resonant mode. Is to be estimated. 
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Rearrange eqn (5-15): 


( 

I 

* 

y 

i 

\ 



On the right hand side of eqn (5-16), the n<ode response power is 
measured, and Q are known, and r is computed by Program RAND. 
So a theoretically correct constant equivalent Input spectrum can 
be identified in simple fashion near the resonant mode. At non- 
resonant regions the problem of input loaded spectra does not 
exist and input identification can proceed in standard fashion 
using the input identification mode of Program RAND. 

5.7 Conclusion to Part II 


Implementation of the RAND procedure in an experimental 
situation brings in a wnole new set of problems: measurement 
accuracy, resolution, input loading effects, noise contamination 
of signals, and statistical averaging. Data show that given 
sufficient care and effort it is possible to obtain useful and 
meaningful response and input Identification estimates. 


The phenomenon of input loading proved to be significant 
experimentally and of definite interest theoretically. The 
simple analytical treatment of this phenomenon reflects the 
simple structures which comprised the laboratory set up. Work 
remains to be done to study the extent to which the present input 
loading models apply to complex field situations, for it is easy 


to speculate that input loading can have considerable practical 
significance. Imagine that a structural component undergoes 
vibration testing and then it is placed in a spacecraft for 
developmental flight testing or undergoes further vibration 
testing at a high structural integration level. In either case. 
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Instrumentation Is provided to characterize the base Input level 
to the component. The results of this study have shown that a 
sharp peak can occur In the Input spectrum due to an Input loading 
effect, and that this peak, being slightly but critically detuned 
to the measured resonant frequency of the component, will have 
only a second order effect on the component resonant response. 

Now consider that the sharp Input peak Is discovered by an envi- 
ronmental specification writer, who routinely draws an envelope 
over U. An erroneous coincidence between the Input peak and the 
component resonance Is thereby guaranteed, and the resulting speci- 
fication Is substantially overestimated. 

Finally, two considerations are added to give prospective to 
the experimental study described In Part II. First, we note that 
all of the experimental work was done without the benefit of a 
modern modal analysis system and without automatic data acquisition 
equipment. The former tool would have greatly improved the accuracy 
of the dynamical modelling of the structure. The latter tool, 
utilized to interface the FFT analyzer output directly to the 
computer, would have greatly Increased the speed and ease of I/O 
operations with Program RAND, thereby making possible the use of 
the full resolution capabilities of the FFT analyzer. The second 
consideration Is simply the recognition that It Is Important to 
distinguish between experimental work done in the laboratory, and 
experimental work done on much more complicated structures under 
field conditions. Obviously success In the laboratory does not 
guarantee success in the field. 


97 


Hopefully the techniques developed and knowledge gained In 
this study, coupled with Improved hardware tools, will add to 
the success with which structural response and Input Identifica' 
tion estimates are made under real field conditions. 
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6.0 Sutnmary 

This study has focused on the problem of the response and 
Input Identification of a structure subjected to combined base 
and acoustic excitation. The study has Included both theoretical 
and experimental approaches. Specific accomplishments are: 

1. The development of a computer program RAMO to make 
response and Input Identification estimates for a 
structure subjected to combined base and acoustic 
excitations. 

2. The development of parametric design guidelines 
for plates subjected to random base or acoustic 
excitation. 

3. Through a combination of analysis and experimenta- 
tion, the evaluation of the applicability of the 
RAND procedure to a laboratory test situation, 

and the identification of the crucial considerations 
which affect the quality of structural Input and 
response estimates. 
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A. 1.1 Case 1: Distant Sound Source 

Refer to Figure 5-10: Model for Input Loading Analysis. Remove 


loudspeaker and acoustic transmission line, and set plate fixture 
displacement = 0. The pressure PQ(t) then acts on the effective 
area A for the first mode of the plate. The equation of motion Is 


y., y„ + k„ y„ = -p.A 
P P P P P P 0 


(A-1) 


which leads to the transfer function 


p p p 


(A-2) 


Next, relate pressure p^ with the Incident pressure p^ emanating 


from the distant sound source: 


P„ = 2 p, + pcy. 


(A-3) 


that Is, the total pressure p^ equals the Incident and reflected 
wave pressures for the blocked plate plus the radiated pressure due 
to the plate velocity y^. Combine eqns (A-2) and (A-3) to obtain 


n /9 n - 

Po/2 P, - 

-w +j2CpiWpU+u)p' 


(A-4) 


where Cp' = pcA)/mp 


A-1 


A. 1.2 Case 2: Near Sound Source 

Refer again to Figure 5-10. The system consists of the loud- 
$peaker« acoustic transmission line, and plate. The plate fixture 
displacement - 0. To facilitate the analysis we redraw the 
system and introduce a displacement function c(x,t) for the acoustic 
transmission line as shown in Figure A-1. 





P. A.C 


”'p' V '‘p 


Figure A-1 Plate Excited by Near Sound Source 

Assume that the dynamics of the air column can be described by 
the one dimensional wave equation 

4 ' 4 (A-5) 

3x^ c^ ar 

subject to the boundary conditions 

4(0, t) = 4,(t) = -y^(t) 

4(L,t) » ^2^t) » -yp(t) 

"’p h ^ ‘’p ^2 ^2 ^ ^ ^ ^ 

A-2 


(A-6) 


For a unit harmonic Input the steady state response at x can be 
written as 

C(x,t) = Rg (H(x.a,)e‘J“^l (A-7) 

The transfer function H(x,u) can be evaluated by substituting 
eqn (A-7) Into eqn (A-5) and then applying the boundary conditions 
(A-6)i 


H(x,w) * cos ux/c + f(u) sin wx/c 
where f(iu) = pAcu.s1nu>L/c-(-u)^ mp+ju)bp+kp)cos uL/c 

T - - 

(-U) m +jwb +k^)s1n wL/c+oAcu cos wL/c 
P P P 


We are Interested In the gage pressure 





x=L 


(A-8) 


(A-9) 


which acts on the plate and the force f which acts on the loudspeaker. 

a 

The loudspeaker force Is obtained by requiring dynamic equilibrium 
for the forces which act on the loudspeaker mass: 

% ■ "’a * *>> h 

where p, = -rc^ 

For the case of harmonic motion 

Pg = pcu (-sin ta)L/c + f(u) cos wL/c) (A-11) 

2 

f, = (-pcAuj f(ui) -w m. + ju) b + k ) 
a a a a 


A- 3 


(A-12) 


The transfer function Mhich relates the pressure Mhich acts on the 
plate to the magnetic force which acts on the loudspeaker Is then 
obtained by dividing (A- 11 ) by (A- 12 ): 


H'h 


•pCu) 


( -sin uiL/c ♦ f(u) cos uL/c j 
-pcwA f(w) - u 


(A- 13 ) 


Substitute the expression for f(u) and simplify: 


P2^^a * 


pcu (-U i"p+jubp+kp) 


(A- 14 ) 


where 




[-(pcAui)^ + (-u^ro.+.lwb +k )(-u^m +jwb +k )l$ 1 n wL/c 

« o a P P P 

+ pcAui [-w^(m^+mp)+ju(b^+bp)+(k^+kp)]cos u»L/c 


A-^ 
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